Kondo effect in quantum dots coupled to ferromagnetic leads with noncollinear 

magnetizations 
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Non-equilibrium Green's function technique has been used to calculate spin-dependent electronic 
transport through a quantum dot in the Kondo regime. The dot is described by the Anderson 
Hamiltonian and is coupled either symmetrically or asymmetrically to ferromagnetic leads, whose 
magnetic moments are noncollinear. It is shown that the splitting of the zero bias Kondo anomaly 
in differential conductance decreases monotonically with increasing angle between magnetizations, 
and for antiparallel configuration it vanishes in the symmetrical case while remains finite in the 
asymmetrical one. 
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Introduction Electronic transport through quantum 
dots attached to nonmagnetic leads was studied exten- 
sively in the past decade - both theoretically and exper- 
imentally. Quantum dots, particularly those including 
only a few or few tens of electrons offer a unique op- 
portunity to study single-electron charging effects and 
electron-electron correlations. In particular, the many- 
body Kondo phenomenon was predicted and then ob- 
served experimentally at low temperatures S H- The 
phenomenon is understood rather well 0, IE IE 0> IE IE 
110|,|11|, although some of its features are still a subject 
of some discussion in the literature - for instance, the 
magnitude of the magnetic-field-induced splitting of the 
Kondo peak in the differential conductance [l2L Ha , fl4| . 

It has been shown recently that the Kondo anomaly 
in electronic transport through a quantum dot at- 
tached symmetrically to ferromagnetic leads becomes 
suppressed in the parallel configuration 0, 0, 0, 0] , 
while in the antiparallel configuration it has features simi- 
lar to those in the case of quantum dots attached to non- 
magnetic leads. This behavior has also been observed 
experimentally [l9j|. The partial or total suppression of 
the anomaly is a consequence of an effective exchange 
field originating from the coupling of the dot to ferro- 
magnetic electrodes. The exchange field gives rise to a 
nonzero spin-splitting of the equilibrium Kondo peak in 
the density-of-states (DOS) in the parallel configuration, 
whereas no such a splitting occurs for antiparallel align- 
ment (in symmetrical situations the overall coupling of 
the dot to the leads is independent of the spin orientation 
for antiparallel alignment). 

The splitting and suppression of the Kondo anomaly 
in electronic transport through quantum dots attached 
to ferromagnetic leads was studied up to now only for 
collinear magnetic configurations (parallel and antipar- 



allel). The noncollinear configuration was considered in 
Ref. [2jj. However, the approximations used to derive 
the appropriate Green functions do not take into account 
the exchange-induced splitting of the anomaly. There- 
fore, in this letter we address the problem of angular 
variation of the zero bias anomaly, i.e., variation of the 
Kondo peak with the angle between magnetic moments 
of the leads. We consider two situations; a symmetrical 
one - like that described above - and a nonsymmetrical 
case corresponding to different strengths and spin asym- 
metries of the coupling to the two electrodes. 

We show that the the Kondo peak in differential con- 
ductance is suppressed already at small deviations from 
the antiparallel configuration. For larger deviations the 
anomaly varies rather slowly with the angle between 
magnetic moments. 

Model and method We consider a single-level quan- 
tum dot described by a Hamiltonian of the general form 
H = H L + H R + H D + H T . The terms Hp describe 
the left (/3 = L) and right ((3 = R) electrodes in the 
non-interacting quasi-particle approximation. The term 
Hjj stands for the dot Hamiltonian and is of the form, 
Hd = J2a £ cr d+dp + Ud^djd^di, where e a is the energy 
of the dot level (we assume that the energy level eo of the 
isolated dot is independent of the spin orientation), and 
U is the Hubbard-type electron correlation parameter. 
The tunneling part, Ht, of the Hamiltonian describes 
electron transitions between the dot and electrodes, and 
is written in the form 1211 
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where cjtg s is the creation operator of an electron with 
the wave vector k and spin s (s — + for majority elec- 
trons and s = — for minority ones) in the lead f3, do- 
annihilates an electron on the dot with spin a [a =T, j 
along the quantization axis of the dot), and is the 

corresponding matrix element. The matrix W^p can be 
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written explicitly as 



Wk/3 



T k/3+ cos(^/2) -T k/3+ sin(^/2) 
n k/3 _ sin(^/2) T k/3 _ oob(^/2) 



(2) 



where 0^ is the angle between the local spin quantiza- 
tion axis in the lead (3 and the quantization axis on the 
dot, and Tkfj s arc the matrix elements for electron tun- 
neling from the dot to the spin majority (s = +) and 
spin minority (s = — ) electron bands in the lead [3 when 
<j)p = 0. As the quantization axis for the dot we assume 
orientation of the exchange field created by the ferromag- 
netic electrodes. For a symmetrical and unbiased system 
one can assume 4>R = — 0l = 9/2, where 8 is the angle 
between magnetic moments of the leads. 

Electric current / has been calculated from a formula 
derived by Meir et al 5] for nonequilibrium situation. 
This formula has been extended to noncollinear magnetic 
configurations. It includes retarded (advanced) and cor- 
relation Green functions of the dot in the presence of bias 
voltage applied to the leads (nonequilibrium situation). 
To calculate the retarded Green functions we used the 
equation of motion method and restricted calculations to 
the large U limit. The higher order Green functions have 
been decoupled according to the procedure described by 
Meir p|. Such an approach gives reasonable results for 
temperatures comparable to the Kondo temperature Tk ■ 
However, we go beyond this approximation by calculating 
the Green functions, renormalized dot level, and occu- 
pation numbers self-consistently. Strictly speaking, the 
bare (spin-degenerate) dot level eo in the self energy Si 
was replaced by the corresponding level spin-split by the 
effective exchange field 22] 



Bexch = (l/ffMs) V n^Re 
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where is a unit vector along the magnetic moment of 
the /3-th electrode, tq is a relevant relaxation time 

is the Fermi distribution functions for the electrode 
0, and It = 2nJ2 k \T kl3s \ 2 S(E - e kf3s ) for « = +,-. 
In the following rf are assumed to be independent of 
energy within the electron band and zero otherwise; 

= r„ (1 + sjo^), where pp is the spin polarization of the 
lead (3 18]. The renormalization of the dot level plays 
an important role in tranport through quantum dots at- 
tached to ferromagnetic electrodes, and leads to the spin 
splitting and suppression of the Kondo anomaly in differ- 
ential conductance 0, fill IrH ITsI ] . The approach based 
on the equation of motion with the renormalization pro- 
cedure including the exchange field gives results which 
are fully consistent with those obtained for collinear con- 
figurations by the perturbational real-time diagrammatic 
approach |23|. 

For U — * oo, the equation of motion for the Green 
function G leads to the Dyson equation (I— go£)G = go, 
where g aa> = 5 aa >(E - e ) _1 and £ = g^ 1 - n^g^ 1 - 
fi _1 S. Here, h„„ = 1 — n,aa and hn-ff — rirrrr, where a = 



—g and n aa i — (d+fV). The self-energy £ is defined as 
S = So + Sr, where IW = ~(*/2) E/3 r L with r n = 
cos 2 (0^/2) + sin 2 (00/2), rfj = sin 2 (0,3/2) + 

cos 2 (0^/2), and = (l/2)(r£ - if) sin^). In 
turn, Si is given by 



de 



2?r E - e n 



ih/rg ' 



2tt E-e 



(4) 



(5) 



The corresponding lesser Green function G < has been 
calculated from the Keldysh equation G < = G 7 "E < G a 
with the self energy calculated from the Ng ansatz. 
Writing S' -S Q = -iT cS and Eg -Eg = -iT = -i(T L + 
F R ), the formula for electric current takes the form 

1 = h I ^ Tr { rLGrf + T R G r T L G a }(f L - f R ), 

(6) 

where T 13 = r^r-ipcS , and e is the electron charge. 

Numerical results Let us write T$ = oFq = aFo. Nu- 
merical calculations have been performed for eo = —0.35, 
kT = 0.001, and To = 0.1 (energy is measured in the 
units of D/50, where D is the width of electron band 
extending from ~D/2 to D/2). Consider first the sym- 
metrical situation: a = 1 and pl = Pr = P- In Fig.l we 
show the corresponding DOS for unbiased (a) and biased 
(b) system for parallel (8 = 0), antiparallel (8 = ir), and 
perpendicular (8 = tt/2) magnetic configurations. In the 
unbiased system, the spin splitting of the Kondo peak 
in DOS is maximal in the parallel configuration and de- 
creases with increasing angle between the magnetic mo- 
ments, vanishing exactly for the antiparallel configura- 
tion. This splitting is a direct consequence of the ex- 
change field exerted on the dot by ferromagnetic elec- 
trodes. The largest exchange field occurs in the paral- 
lel configuration, so the corresponding splitting of the 
anomaly is also maximal in this configuration. In the an- 
tiparallel configuration the exchange fields from the two 
electrodes cancel each other so the splitting vanishes (see 
Fig. 1(a)). A further consequence of the spin splitting of 
DOS is a suppression of the zero bias anomaly in the 
differential conductance, as will be discussed below. 

In the biased system, the Kondo peak in DOS be- 
comes split due to different electrochemical potentials of 
the two electrodes. Such a splitting exists also in non- 
magnetic situations. The exchange field in parallel and 
noncollinear configurations introduces further splitting, 
as shown in Fig. 1(b), similarly as an external magnetic 
field does in the case of quantum dots attached to non- 
magnetic leads |l2^ . 

Numerical results for the corresponding differential 
conductance and tunnel magnetoresistance (TMR) vs 
bias voltage are shown in Fig. 2 for several values of the 



3 



angle 9 between magnetic moments of the leads, and for 
two values of the spin polarization factor p. In the an- 
tiparallel configuration (6 = it) there is a well defined 
zero bias anomaly due to the Kondo correlations. When 
the angle between magnetic moments decreases, the con- 
ductance peak becomes split and its intensity becomes 
suppressed. The splitting as well as the suppression of the 
intensity grow with decreasing 9. Finally, in the parallel 
configuration only low-intensity peaks occur at nonzero 
positive and negative bias. The results are consistent 
with the limiting situations of collinear c onfigura tions, 
9 = and 9 = tt, considered in Refs [H |H Elll . It is 
interesting to note a fast decrease of the peak intensity 
already at small deviations from the antiparallel configu- 
ration, as follows from Fig. 2(a) and (b). Variation of the 
peak intensity with 9 is rather small for 9 < ir/2. Com- 
parison of Fig. 2 (a) and (b) also shows that the splitting 
increases with increasing spin polarization of the leads, 
while the peak intensity decreases with increasing p. 



at V — 0. Two side local minima of TMR occur then 
at some nonzero voltages. The maximum at V = in- 
creases and becomes more flat when the angle between 
magnetic moments decreases from 9 = tt to 9 = 0. 
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FIG. 1: Density of states in the Kondo regime in (a,b) sym- 
metrical (a = 1, p = 0.2) and (c,d) asymmetrical (a = 0.1, 
Pl = 0.4, pa = 0.8) systems for three indicated values of the 
angle between magnetic moments and for unbiased (a,c) and 
biased (b,d) systems. The other parameters as described in 
the text 

The associated TMR ratio, defined by the system re- 
sistance R{9) as [R(9) - R{9 = 0)}/R(9 = 0), shows a 
peculiar behavior displayed in Fig. 2 (c) and (d) for two 
values of p. For 9 = tt the TMR has a dip in the zero 
bias regime and is negative. The negative sign is a conse- 
quence of the suppression of the Kondo peak in the par- 
allel configuration, which makes the system more con- 
ducting in the antiparallel configuration (inverse TMR 
effect). The dip at zero bias disappears when the initial 
configuration departs from the antiparallel one, and in- 
stead of a minimum there is a local maximum of TMR 
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FIG. 2: Bias dependence of the Kondo anomaly in differential 
conductance (a,b) and of the TMR effect (c,d), calculated for 
indicated values of the angle 9 between magnetic moments 
of the leads and for two different polarization factors p. The 
other parameters as in Fig.l 

In the situation described above the quantum dot was 
coupled symmetrically (in the parallel configuration) to 
the leads. When the coupling is not symmetrical, then a 
certain splitting and suppression of the zero bias anomaly 
exist also in the antiparallel configuration [2^, 0] ■ This 
is because the exchange fields from the two leads do not 
compensate each other as they did when the dot was 
coupled symmetrically. Accordingly, the splitting of the 
Kondo anomaly may be observed also when only one of 
the leads is ferromagnetic while the other one is nonmag- 
netic 0]. 

Let us consider now angular variation of the Kondo 
anomaly in a strongly asymmetric system; a = 0.1, 
Pl = 0.4, pb = 0.8. The corresponding DOS is shown in 
Fig. 1 (c) and (d) respectively for the unbiased and biased 
system, and for parallel, antiparallel and perpendicular 
magnetic configurations. At V = the Kondo peak in 
DOS is spin-split in the antiparallel configuration and the 
splitting increases when the configuration varies from an- 
tiparallel to parallel. In the biased system the peaks are 
additionally split due to different electrochemical poten- 
tials of the two leads, but intensity of one of the two com- 
ponents of the peak is relatively small and only weakly 
resolved in Fig. 1(d). The asymmetry in peak intensities 
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FIG. 3: Bias dependence of the Kondo anomaly in differential 
conductance (a) and TMR (b), calculated for indicated values 
of the angle and lead polarization factors. 



is a consequence of the spin asymmetry in the coupling 
of the dot to metallic electrodes, and a difference in spin 
polarizations, pl ^ Pr- The differential conductance and 
TMR are shown in Fig. 3 (a) and (b), respectively. In the 
parallel configuration the dominant Kondo anomaly oc- 
curs in the spin-up channel and for eV > (electrons flow 
from left to right). The corresponding peak at eV < 



is significantly weaker. In the antiparallel configuration, 
on the other hand, the Kondo peak in differential con- 
ductance occurs for eV < only. This behavior can 
be accounted for by taking into account behavior of the 
Kondo peaks in DOS and spin asymmetry in the coupling 
to electrodes pij . 

The corresponding TMR is shown in Fig.3(b). It is 
interesting to note that TMR is highly asymmetric with 
respect to the bias reversal. It becomes positive for eV 
exceeding a certain positive value, and negative below 
this voltage. This is a consequence of the fact that for 
positive eV the Kondo peak in differential conductance 
is clearly visible in the parallel configuration, whereas for 
eV < the Kondo peak occurs in the antiparallel config- 
uration. Such a behavior of the conductance and TMR 
may be interesting from the point of view of applications 
in mesoscopic devices like diodes, for instance. 

In summary, we have analyzed the Kondo phenomenon 
in a quantum dot attached to ferromagnetic leads with 
noncollinear magnetizations. Our results show a mono- 
tonic dependence of the Kondo peak suppression with 
the angle between magnetic moments of the leads. A 
fast decrease of the Kondo peak takes place already at 
small deviations from the antiparallel configuration. 

This work was supported by EU trough RTNNANO 
contract No MRTN-CT-2003-504574, RTN 'Spintronics' 
HPRN-CT 2002-00302. 
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